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Abstract. A fermion model consisting of two chains with interchain tunneling is 
formulated and solved exactly by the Bethe ansatz method. The interchain tunneling 
leads to Cooper pair like bound states and a threshold energy is required to overcome 
the binding energy. The correlation functions manifest superconducting properties. 
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Since several years there has been considerable interest in unusual mechanisms 
for the formation of Cooper pairs because of the possible relevance for high T c 
superconductors. A simple model showing a pairing transition was recently proposed 
in [1] to explain the properties of the high-temperature cuprate superconductors. 
The underlying mechanism leading to the superconducting transition is an interlayer 
tunneling phenomenon. Although the model in [1] is two-dimensional, it is quite 
interesting to investigate its one-dimensional variant, as it has been conjectured [2] 
that properties of the one- and two-dimensional models have common aspects. Exact 
results in one dimension are often more easily available than in two dimensions and 
may provide a testing ground for approaches intended for more complex problems. 

In this paper we present the exact solution to a new integrable model consisting of 
two chains with an interchain tunneling which may be considered as the simplest one- 
dimensional variant of the model in [1]. It will be shown that the attractive interaction 
leads to bound states of Cooper type pairs and the model has superconducting 
properties. 

The Hamiltonian under consideration is the following 

L 

TC = -^2v[c+ lCj+hl + c+ 2 c j+h2 - Vc+ lC + +hl {cj2C j+1 , 2 + Cj-i, 2 c j2 ) + h.c] V, (1) 

where the Fermi operator Cj a creates an electron at site j with a sublattice index 
a(a = 1,2), L is the number of lattice sites, and h.c. denotes the Hermitian 
conjugate. We impose the restriction of no simultaneous occupancies of sites on 
different sublattices which are nearest neighbours, i.e. (j, 1); (j, 2) and (j, 1); (j — 1,2), 
and V is the projector on the subspace of allowed states. The Hamiltonian (1) contains 
two different terms. The first describes the hopping of individual electrons with 
transition rate set equal to 1, and the second describes the hopping of pairs of electrons 
with rate V (> 0). The competition of these terms is the main property of the model 
under consideration. It should be noted that this is similar to the one-dimensional 
model proposed by Penson and Kolb [3]. We hope that the study of (1) will also help 
to understand this nonintegrable model. 

The exact solution for the eigenstates and eigenvalues of Hamiltonian (1) can be 
obtained within the framework of the Bethe ansatz method [4-6]. The structure of 
the Bethe ansatz equations follows from the solution of the two-particle problem. The 
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nonvanishing elements of the S'-matrix are 



Sim = = — ^re- tt , ( 2 ) 

sm(k + %r\) 

sim=s\m=^ i \ T/=-inv. 

A necessary and sufficient condition for the applicability of the Bethe ansatz method 
are the Yang-Baxter equations [4]. Since the scattering matrix (2) satisfies these 
equations, our model (1) is integrable. 

The discrete Bethe ansatz equations are derived following the standard procedure 
[7] by imposing periodic boundary conditions. Each state of the Hamiltonian is 
specified by one set of charge rapidities kj representing the momenta of the electrons 
and one set of additional rapidities A a . All rapidities within a given set have to be 
different, corresponding to Fermi statistics. These rapidities are determined by the 
Bethe ansatz equations 



ik-L TT \ i n i \1 TT sin [U k 3 ~ A «) + § V] 

e*^ = exp -(ki-h) — A--^ ; f— 



1=1 

N ■ rl / a 7 \ , r l M 



(3) 



-j-j- sin[|(A a - kj) + %rf[ yr sin[^(A a - A /3 )+ C 77] 
ji sin[i(A a - kj) - §t/] ~ Al sin[|(A a - A^)- C 77] ' 

The energy and the momentum of the corresponding state are given by 

N N 

£=-2^>os% P = J2 k J- ( 4 ) 
j'=i i=i 

We remark that in contrast to [4-6] in our model the electron number on each 
separate sublattice is not conserved. In the case under consideration the conserved 
quantities are the total number of electrons iV and the number M of pairs of 
consecutive electrons which are on the same sublattice. 

Each eigenstate of the Hamiltonian is specified by sets of rapidities which have 
to satisfy the Bethe ansatz equations (3). The structure of the solutions of these 
equations is different for 77 < and 1] > 0. In the case 1] > all kj are real and 
particles move independently. There are no Cooper pairs in the system. Equations 



(3) are reduced to one set of equations for k 



j 

N N 

2 



kjL - - J2( k 3 ~ k i) + J2 ~ ki "> = 27vI J> 



i=i i=i 

oo 



(5) 

§(k) = - + V ex P(-^) sm(uk) 
1 9 ^ z/ cosh(^)' 



i/=i 



where /j are integer (half-integer) numbers for odd (even) N. (We have included the 
linear term k/2 in the definition of the function $(/c) such that it is identical to the 
phase function used in [10].) In the thermodynamic limit eqs. (4,5) are replaced by 
the integral equation for the distribution function of particles p(k) 

2np(k) - [ ° <p(k - k')p{k')dk' = 1 - -p, <p(k) = &(k) (6) 
J-k 2 

with subsidiary condition for the particle density p 

-k 



p{k)dk = p = N/L. (7) 

— k 

In the opposite case rj < (V > 1) pairs of electrons are energetically more 
favourable as compared to unpaired particles. In the ground state we only have Cooper 
pair like bound states with finite pairing energy which are characterized by pairs of 
complex wavenumbers 

kt = A a ±c\ V \. (8) 

Inserting (9) into the Bethe ansatz equations (3) we obtain the following set of 
equations 

M M 

2LA a - 2 J^(A Q -A p )-J2 ©( A a " A /3! M) = J a , 

(3=1 (3=1 

" (9) 
E 1 = ^eo(A Q ), eo(A) = — 4cosh?]cosA, 

a=l 

B(A;r]) = 2 arctan(coth?ytan A/2), —n < ©(A;??) < 7r, 
where J a are integers or half-odd integers depending on the parity of M . In the 
thermodynamic limit eqs. (10) lead to the integral equation for the distribution 
function of Cooper pairs a (A): 

2yra(A) + 6' (A - A'; \ V \)a(A')dA' = 2 (l - ip) , (10) 



/An 
-An 



with subsidiary condition 

a(K)dK=\p. (11) 

A 1 

In order to understand the superconducting properties of the model, we shall 
investigate the long-distance behaviour of the correlation functions. For this purpose 
we shall use the methods of two-dimensional conformal field theory [8,9]. Our model 
has just one type of gap less excitations which correspond to "holes" and "particles" 
in the /-distribution or J-distribution for the case rj > and rj < 0, respectively. 
Therefore, an application of this theory presents no difficulties and can be carried 
through as in [10,11]. The results of these calculations are the following. 

For rj > the long-distance power-law behaviour of the density-density correlation 
functions is given by 

(n(r)n(O)) ~ p 2 + Ar~ a cos(7rpr), n(r) = ^ c+ T c rr . (12) 

T = l,2 

The critical exponent a is expressed in terms of the dressed charge Z which in our 
case is given by 

a = Z 2 /2, Z = 2irp(k ). (13) 

The superconducting properties of the system manifest themselves in the behaviour 
of the pair correlation functions for which we obtain 
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G p (r, r) = (c+c+ +1 lCOr c lr ) ~ Br-**, f3 = 2/Z 2 = -. (14) 

a 



In the same way we can consider the case rj < 0. The expressions (13)-(15) also 
hold in this case, but the dressed charge Z is given in terms of the distribution function 

Z = 27rt(A ). (15) 

Comparison of eqs. (13)-(16) with the corresponding results for the model of 
electrons with correlated hopping [11] permits us to obtain a simple relation between 
the critical exponents of these models 

P(P, V)=(l~ \f>} /3(c) V ) ' V > °' 

p(p,ti)=(i-\p) p(c) (j^TpAviy v<o, (l6) 

< p < 1, 



where (3^ is the critical exponent of the superconducting correlation function for 
electrons with correlated hopping [11]. Numerical results for the critical exponent (3 
are shown in Fig. 1. 

In addition to (17) there is a relation of the charge stiffnesses D(p, rf) and _D( C ) (p, rf) 
of both models 



D{p,r 1 )=[l-\p^e^D {c) (^-^,r 1 ) : r, > 0, 



D(p,r,)=(l-±p)D (e) (^,\v\), V<0, (17) 
< p < 1. 

Figs. 2 and 3 show numerical results for the charge stiffness and transport mass of 
the charge carriers defined by 

m e D 

where m e is the bare mass and D° = ^ sin ^p is the charge stiffness of free fermions. 
For densities p approaching the maximum allowed value 1 the masses diverge due to 
a (trivial) metal-insulator transition. 

The existence of a simple relation between the critical exponents, a/3 = 1, is 
certainly a consequence of the universal Luttinger behaviour of critical models with 
one type of gapless excitations [12]. As our model of interchain tunneling (1) is 
described by a conformal field theory with central charge c = 1, the critical exponents 
depend in a nonuniversal way on the particle density p and the interaction strength r\. 
Although the model of electrons with correlated hopping [11] is different, there exists 
a mapping of the phase diagram of the first model to the latter one with a simple 
relation of the phase points and the corresponding values of the exponents (3 and /?( c ) , 
respectively. Note that the mapping is one-to-one, however, the physical properties of 
the models differ. In contrast to [11] we see that for (1) with couplings r\ > there 
are no densities p for which (3 < a. This means that the dynamics is dominated by 
the motion of individual electrons rather than by the motion of pairs, and we have no 
"superconducting" phase for any particle concentration. 

For rj < the situation is quite different. As in the model with correlated hopping 
there is a regime for which (3 < a, see. Fig. 4. In this case the pair correlation function 
has a slower decay than the density-density correlation function and thus dominates. 
In one-dimensional systems which do not undergo a condensation and do not have 



long-range order such behaviour of correlation functions is the closest analogy to the 
existence of a superconducting phase [13]. 
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Figure captions 

Figure 1. (a) Depiction of the critical exponent (3 in dependence on the particle 
density p for different values of r\ =0.1, 0.5, 1, 10. (b) The same as for (a) with 
negative values of r\, \r}\ = 0.05, 0.1, 0.5, 1, 10. 

Figure 2. Charge stiffness D c and effective mass m of the charge carriers for different 
interaction strengths r\ = 0.1, 0.5, 1, 2, 10. 

Figure 3. Charge stiffness D c and effective mass m for negative values of \rj\ = 
0.1, 0.5, 1, 2. 

Figure 4. Depiction of the phase diagram for r\ < 0. The solid line separates regions 
with dominating pair and dominating density-density correlations ((3 < a and j3 > a), 
respectively. For large r\ the critical density approaches (4 — 2y/2)/3 = 0.39...). 
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